Thermodynamics of the one-dimensional s=l XXZ Heisenberg modeh 

analytical results 
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We apply the results derived by Rojas et a/. |l| to derive the /3-expansion of the Helmholtz free 
energy of the spin-1 XXZ Heisenberg model up to order The analytical expansion obtained is 
valid for all phases of this model. Our curves of the specific heat fit well Blote's numerical resultsQ 
in the high temperature regime. 
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The spin-1 XXZ Heisenberg model is non-integrable 
and therefore cannot be solved by the algebraic Bethe 
ansatz method Since the 70 's its thermodynamics has 
been studied numerically |0, 4. 0. More recently, 

Yamamoto and Miyashitala, 6, u\ applied a Monte Carlo 
method to study the specific heat, static magnetic suscep- 
tibility and magnetization of rings and chains of different 
sizes as a function of temperature and extended those re- 
sults to the thermodynamic limit. Bao et al.^ derived 
a set of self-consistent equations to approximately ob- 
tain the thermodynamics of this model for all range of 
temperature, and solved those equations numerically. 

The cummulant series |P yields a /3-expansion of the 
Helmholtz free energy for exact integrable as well as 
for non-integrable models. Recently Rojas et al.^ 
showed for any one-dimensional chain model with peri- 
odic boundary condition, invariance under spatial trans- 
lation and interaction between nearest neighbours that, 
in the thermodynamic limit, the coefficient of the high 
temperature expansion of arbitrary order /3" can be de- 
rived from an auxiliary function {p. 

Our aim in this work is to apply the results of Ref. |l| 
to derive an analytical /3-expansion of the Helmholtz free 
energy per site of the anisotropic spin-1 XXZ Heisenberg 
model with single-ion anisotropy, up to order /3^. 

The hamiltonian of the spin-1 XXZ Heisenberg model 
with anisotropy pQl is 
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where N is the number of sites in the periodic chain, 
A the anisotropy constant, h the external magnetic field 
in the z-axis and D the single-ion anisotropy parameter. 
We have 5"=^ = ^{Sf ±i S^). do not fix the sign of 



any constant in eq.(|I|). 

Applying the results of Ref. |^ to the spin-1 XXZ 
Heisenberg model with single-ion anisotropy, we obtain 
its Helmholtz free energy W(/3) up to order (3^: 
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We point out that this analytic expansion is vahd for any 
arbitrary set of values of the parameters J, A, D and 
h in the high temperature regime. Therefore it applies 
equally well to any phase of the one-dimensional spin-1 
XX Z Heisenberg model||ll|, [ij, |l3l. The coefficients for 
each order in f] in the expansion (Id) are exact. 

In addition, from eq.(^ we see that the thermody- 
namic properties of this model are insensitive to the sign 
of the constants J and h. From now on we take J > 0. 
By redefining all the parameters of the hamiltonian (|^), 
eq.(|^) becomes an expansion in powers of (J/?)". 

A beautiful numerical work by B16te[^ tabulates nu- 
merical values for the specific heat at varied temperatures 
(including the high temperature region) for the s = 1/2 
and s = 1 Heisenberg models. In all his calculations 
Blote has set h = 0. To verify the correctness of our ex- 
pansion (|^) we compare the analytical expression for the 
specific heat derived from it with the numerical results of 
Ref. ^. Comparison of the hamiltonian (1) in Ref. || and 
ours, given by eq.(|l|), yields the following correspondence 
of parameters: Jy = —A/2 and J± — — J/2. In Fig. |l|we 
compare our analytical expression for the specific heat 
derived from eq. for anisotropic Heisenberg models in 
the absence of single-ion anisotropy {D = 0) with dif- 
ferent values of J in the antiferromagnetic phase. In the 
high temperature region our curves fit well the numerical 
data. 

It is certainly important to verify that eq.(|^) applies 
equally well to the ferromagnetic phase of the spin-1 
Heisenberg model. To show this we have Fig. |^ where 
we plot the specific heat, in the high temperature region, 
of the isotropic Heisenberg model with different values for 
the single-ion parameter D in this phase. As in Fig. 
our curves in Fig. |2| fit pretty well Blote's numerical re- 
sults in the high temperature region. 

Our analytical results for the specific heat, obtainable 



from (g), fully agrees in its isotropic limit with that of 
section 5.2 of Ref. || in the high temperature limit. We 
also recover the high temperature static magnetic sus- 
ceptibility for the isotropic Heisenberg model [eq.(4.5a) of 
Ref. 1^ . However, we do not agree with the high tempera- 
ture limit of the correlation function for the z-component 
of spin between nearest neighbors (S'pS'f ) of Ref. |[ The 



correct Hmit is {SgSf) 



A/3. 



In summary, in this letter we applied the method of 
Ref.0 to a non-integrable spin-1 XX Z Heisenberg model 
with single-ion anisotropy in the presence of an external 
magnetic field in the z-direction. We obtained the ana- 
lytical /3-expansion of the Helmholtz free energy per site 
of the model up to order . Each coefficient in expansion 
(^) is exact and valid for any phase of the model. 

Our curves in the high temperature region fit well 
Blote's numerical results for the spin-1 XX Z Heisenberg 
model. We recover the few analytical results known in 
the literature H, |^ about the isotropic spin-1 Heisenberg 
model. We correct the high temperature limit of (SqSi) 
of Ref. ^. For the first time in the literature we have an- 
alytical results easily handled for the anisotropic spin-1 
Heisenberg model with single-ion anisotropy and in the 
presence of an external magnetic field. 

Finally we should mention that our calculation pro- 
cedures have been implemented in the symbolic compu- 
tational language Maple|l^. Currently, refinements are 
being made so that the expansion of the Helmholtz free 
energy (H) for the s = 1 Heisenberg model can be ex- 
tended to higher orders in /3. 

An enlarged version of this letter will be published else- 
where. 
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FIG. 1: Comparison of our analytic expression for tfie spe- 
cific heat and Blote's numerical result [H] for the anisotropic 
Heisenber model with different anisotropies and h = 0. Solid 
lines stand for the specific heat of our analytical result and 
dotted lines correspond to data from Ref. 0. 
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FIG. 2: Comparison of our analytic expression for the specific 
heat and Blote's numerical results|^] for the isotropic Heisen- 
berg model with different single-ion parameters and h — 0. 
Solid lines stand for our results and dotted lines correspond 
to data from [01. 



